Mon. Not. R. Astron. Soc. 414. [2528lf2540l (20111 Printed 22 June 2011 (MN I*TeX style file v2.2) 



The LBA Calibrator Survey of southern compact 
extragalactic radio sources — LCSl 



o 

m 
(N 

o 

u 

6 

(N 
> 

O 

\o 

(N 



Leonid Petrov^*, Chris Phillips^, Alessandra Bertarini^, Tara Murphy^'^, and 
Elaine M. Sadler^ 

^ADNET Systems, Inc. /NASA GSFC, Code 610.2, Greenbelt, MD 20771 USA 

^ CSIRO Astronomy and Space Science, PO Box 76, Epping, NSW 1710, Australia 

^Institute of Geodesy and Geoinformation, University of Bonn, Nussallee 17, Bonn, Germany and 

Max Planck Institute for Radioastronomy, Bonn, Germany 

'''Sydney Institute for Astronomy, School of Physics, The University of Sydney, NSW 2006, Australia 
^School of Information Technologies, The University of Sydney, NSW 2006, Australia 



22 June 2011 



ABSTRACT 

We present a catalogue of accurate positions and correlated flux densities for 410 
flat-spectrum, compact extragalactic radio sources previously detected in the AT20G 
survey. The catalogue spans the declination range [—90°,— 40°] and was constructed 
from four 24-hour VLBI observing sessions with the Australian Long Baseline Array 
at 8.3 GHz. The VLBI detection rate in these experiments is 97%, the median un- 
certainty of the source positions is 2.6 mas, and the median correlated flux density 
on projected baselines longer than 1000 km is 0.14 Jy. The goals of this work are 
1) to provide a pool of southern sources with positions accurate to a few milliarc- 
sec, which can be used for phase referencing observations, geodetic VLBI and space 
navigation; 2) to extend the complete flux-limited sample of compact extragalactic 
sources to the southern hemisphere; and 3) to investigate the parsec-scale properties 
of high-frequency selected sources from the AT20G survey. As a result of this VLBI 
campaign, the number of compact radio sources south of declination —40° which have 
measured VLBI correlated flux densities and positions known to milliarcsec accuracy 
has increased by a factor of 3.5. The catalogue and supporting material is available at 
[http : / / astrogeo . org/lcsl . 
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1 INTRODUCTION 

Catalogues of positions of compact extragalactic radio 
sources with the highest accuracy are important for many 
applications. These include imaging faint radio sources in 
the phase referencing mode, differential astrometry, space 
geodesy, and space navigation. The method of VL B I first 
proposed by iMatveenko. Kardashev fc Sholomitskiil (| 19651 ) 
allows us to derive the position of sources with nanoradian 
precision (1 nrad « 0.2 mas). The first catalogue of source 
coordinates determined with VLBI contained 35 objects 
ICohen fc Shafi^erl Il97lh . Since then hundreds of sources 
have been observed under geodesy and astrometry VLBI 
observing programs at 8.6 and 2.3 GHz (X and S bands) us- 
ing the Mark3 recording system at the International VLBI 
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Service for Geodesy and Astrometry (IVS) network. Analy- 
sis of these observations res ulted in the ICRF catalogue of 
608 sources (|Ma et aLlllOQ^ ). 

The Very Long Baseline Array (VLBA) was later 
used to measure the positions of 4000 compact ra- 



jBeaslev et all 20021: Fomalont et al. 2003 


: Petrov et al.l 


I2OO5I. 20061: Kovalev et al.l 


20071: Petrov et al. 


2007a) and the 


geodetic program RDV ( 


Petrov et all l2009bl). All sources 



with declinations above —45° detected using Mark3/Mark4 
under IVS programs were re-observed with the VLBA in the 
VCS and RDV programs, which significantly improved the 
accuracy of their positions. As a result of these efforts, the 
probability of finding a calibrator with a VLBI-determined 
position greatly increased. In the declination range 5 > —30° 
the probability of finding a calibrator within a 3° radius of 
a given position reached 97% by 2008. 
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Since the VLBA is located in the northern hemisphere, 
observations in the declination zone [—50°,— 30°] are diffi- 
cult and the array cannot observe sources with 5 < —52°. In 
2008, the probability of finding a calibrator within a radius 
of 3° was 75% in the declination zone [-40°, -30°] and 42% 
for declinations south of —40°. The VLBI calibrator lislQ 
in 2008 had 524 sources in the zone [+52°, +90°], but only 
98 objects in the zone [—52°, —90°] which cannot be reached 
by the VLBA. These southern sources were observed dur- 
ing geodetic experiments and during two dedicated s o uther n 
hemisphere astrometry campaigns |Fev et al ] l2004l . l2006h . 
The reason for this disparity is the scarcity of VLBI anten- 
nas in the southern hemisphere, particularly stations with 
dual frequency S/X receivers and geodetic recording sys- 
tems. Also until recently there has been a lack of good all-sky 
catalogues suitable for finding candidate VLBI calibrators. 

The Australian Long Baseline Array (LBA) consists of 
six antennas located in Australia with the South Africa sta- 
tion HARTRAO often joining in. This VLBI network operates 
3-4 observing sessions a year, each about one week long. 
Although the hardware used by the LBA was not designed 
for geodesy and absolute astrometry observations, it was 
demonstrated by IPetrov et al] (|2009al ') that despite signifi- 
cant technical challenges, absolute astrometry VLBI obser- 
vations with the LBA network is feasible. In a pilot exper- 
iment in June 2007, the positions of participating stations 
were determined with accuracies 3-30 mm, and positions of 
five new sources were determined with accuracies 2-5 mas. 

Inspired by these astrometric results, we launched the 
X-band LBA Calibrator Survey observing campaign (LCS), 
with the aim of determining milliarcsecond positions and 
correlated flux densities for one thousand compact extra- 
galactic radio sources at declinations south of —30°. The 
overall objective of this campaign is to match the density of 
calibrators in the northern hemisphere and so eliminate the 
disparity. 

We have three long-term goals in this campaign. Firstly, 
setting up a dense grid of calibrators with precisely known 
positions within several degrees of any target will make make 
phase-refer encing o b servat ions of weak targets feasible. Ac- 
cording to IWrobell (|2009l '). 63% of VLBA observations in 
2003-2008 were made in the phase referencing mode. A 
dense grid of calibrators also makes it possible to do dif- 
ferential astrometry of Galactic objects such as pulsars and 
masers, and allows direct determ ination of the par allax at 
distances up to several kiloparsec (|Deller et al.ll200"9l ). These 
sources form the pools of targets for observations under the 
geodesy programs and for space navigation. 

Our second goal is to extend the complete flux-limited 
sample of compact extragalactic radio sources (with emis- 
sion from milliarcsecond-size regions) to the entire sky. Ac- 
cording to Kovalev (private communication, 2010), analysis 
of the log A*'-log S diagram of the VLBI calibrator list sug- 
gests the sample of radio-loud Active Galactic Nuclei (AGN) 
is complete at the level of correlated flux density 200 mjy 
at X-band at spatial frequencies 25 MX at 5 > —30°. Ex- 
tending this complete sample to the entire sky will make it 
possible to generalize the properties of the sample, such as 
distribution of compactness, distribution of brightness tem- 

^ Available at I http : //astrogeo ■ org/rf c | 
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perature, bulk motion, viewing angle, irregularities of the 
spatial distribution, to the entire population of radio loud 
AGN. 

The third goal is to investigate the properties of 
high-frequency selected radio sources from the Australia 
Telescope 20 G Hz (AT20G) survey (|Murphv et all l2010l : 
iMassardi et al.ll2010l ). Obtaining observations of a subsam- 
ple of AT20G sources with milliarcsecond resolution will al- 
low us to investigate the properties, such as spectral index, 
polarization fraction and variability, of a population of ex- 
tremely compact sources. 

In this paper, we present the results from the first four 
24-hour experiments observed in 2008-2009. The selection of 
candidate sources from the AT20G catalogue is discussed in 
section[51 The station setups during the observing sessions is 
described in section [3] The correlation and post-correlation 
analysis, which is rather different from ordinary VLBI exper- 
iments, is discussed in sections U and [S] An error analysis of 
single-band observations, including evaluation of ionosphere- 
driven systematic errors, is given in subsection l5.2l The cat- 
alogue of source positions and correlated flux densities is 
presented in section [6] and the results are summarized in 
section [T] 



2 CANDIDATE SOURCE SELECTION 

The Australia Telescope 20 GHz (AT20G) survey is a blind 
radio survey carried out at 20 GHz with the Australia 
Telescope Cornpact Array (ATCA) between 2004 and 2008 
iMurphv et al.ll2010h . It covers the whole sky south of dec- 
lination 0°. The source catalogue is an order of magni- 
tude larger than previous catalogues of high-frequency radio 
sources, with 5890 sources above a 20 GHz flux-density limit 
of 40 mJy. All AT20G sources have total intensity and po- 
larization measured at 20 GHz, and most sources south of 
declination —15° also have near-simultaneous flux-density 
measurements at 4.8 and 8.6 GHz. A total of 1559 sources 
were detected in polarised total intensity at one or more of 
the three frequencies. There are also optical identiflcations 
and redshifts for a signiflcant fraction of the catalogue. 

This high-frequency catalogue provides a good starting 
point for s electing bright, compa ct sources and candidate 
calibrators. IMassardi et al.l (|20ld ) show that almost all the 
flat-spectrum AT20G sources with spectral index a > —0.5 
are unresolved on scales of 0.1-0.2 arcsec in size at 20 GHz 
(see their fig. 3.4). The few exceptions are either (i) fore- 
ground Galactic and LMC sources such as planetary neb- 
ulae, HII regions and pulsar wind nebulae, which have a 
flat radio spectrum due to their thermal e mission but are 
usually resolved on scales of a few arcsec (|Murphv et al.l 
l2010h . or (ii) flat-spectrum extragalactic sources which are 
gravitationally-lensed, like PKS 1830-211. 

For 64% of sources from the AT20G catalogue, flux den- 
sities were determined in three bands, 5.0, 8.4 and 20 GHz. 
We used these measurements to calculate the spectral index 
a {F oc f"). We selected a set of 684 objects, not previously 
observed with VLBI, which had a flux density > 150 mJy at 
8.3 GHz and a spectral index > —0.6. We then split this sam- 
ple into two subsets: 410 high priority sources with flux den- 
sities > 200 mJy and spectral indices > —0.5, and all others. 
In addition, we selected 14 flat-spectrum objects that did 
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Table 1. The LBA network stations. The average System Equiv- 
alent Flux Density (SEFD) at 8.3 GHz at elevation angles > 60° 
achieved in LCSl experiments is shown in the last column. 



Name 




<i>gc 




A 


Diam 


SEFD 


ATCA-104 


-30° 


.15 


149° 


.57 


22 


m 


350 


Jy 


CEDUNA 


-31° 


.70 


133° 


.81 


32 


m 


490 


Jy 


DSS45 


-35° 


.22 


148° 


.98 


34 


m 


120 


Jy 


HOBART26 


-42° 


.62 


147° 


.44 


26 
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620 


Jy 


MOPRA 


-31° 


.10 


149° 


.10 


22 
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Jy 
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Jy 
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Figure 1. The LBA network used for LCSl observations. 



not have flux density measurements at 5 GHz and 8 GHz in 
AT20G. Their spectral indices were determined by analyzing 
historical single dish observations found in the Astrophysica l 
CATalogs support System CATS (|Verkhodanov et al.lll997l ') 
database, which by 2010 included data from 395 catalogues 
from radioastronomy surveys. 

In addition to these target sources, we identified a set 
of 195 sources that we used for calibrator selection. These 
were bright sources previously observed at the VLBA or IVS 
network with position known with accuracies better 0.5 mas 
and with correlated flux densities > 500 mjy on baselines 
longer than 5000 km. 



2.1 The Long Baseline Array 

The LBA network used for LGSl consists of 6 stations: 
ATCA-104, CEDUNA, DSS45, HOBART26, MOPRA and PARKES 
(Fig. m Table [!}. The maximum baseline length of the net- 
work is 1703 km, the maximum equatorial baseline projec- 
tion is 1501 km, and the maximum polar baseline projection 
is 963 km. 



2.2 Observation scheduling 

Observation schedules were prepared with the software pro- 
gram sur_sked. The scheduling goal was to observe each tar- 
get source with all antennas of the array in three scans of 
120 seconds each in the first three experiments and in four 
scans in the last experiment. After tracking a target source 
for 120 seconds, each antenna immediately slewed to the 
next object. The minimum time interval between consecu- 
tive scans of the same source was 3 hours. The scheduling 
software computed the elevation and azimuth of each candi- 
date source, with constraints on the horizon mask and max- 
imum elevations to which antennas can point. It calculated 
the slewing time taking into account the antenna slewing 
rate, slewing acceleration and cable wrap constraints. The 
sources were scheduled to minimize slewing time and fit the 
scheduling constraints. A score was assigned for each tar- 
get source visible at a given time, depending on the slewing 
time, the history of past observations in the experiment and 
the amount of time that the source would remain visible to 
all antennas in the array. 

Every hour a set of four calibrator sources was observed: 
two scans where all antennas had an elevation in the range 
[12°, 45°], one scan at elevations [32°, 45°] and one scan at 
elevations [45°, 90°]. Scans with calibrator sources were 70 
seconds long. The scheduling algorithm for each set found all 
combinations of calibrator sources that fell in the right ele- 
vation range, and selected the sequence of four objects that 
minimized the slewing time. The purpose of these observa- 
tions was 1) to serve as amplitude and bandpass calibrators; 
2) to improve the robustness of estimates of the path delay 
in the neutral atmosphere; and 3) to tie the positions of new 
sources to exist ing catalogues such as the ICRF catalogue 
(|Ma et al.ll 19981 ). 

The scheduling algorithm assigned four scans for 18% 
of target sources, three scans for 76% of sources, two scans 
for 3% of sources, and one scan for 3% of sources. Of the 26 
targets with fewer than three scans in a given observing ses- 
sion, 18 were observed in two experiments. During a 24 hour 
experiment, 11.5-12.0 hours were used for observing target 
sources, 1.5-2.0 hours for observing calibrators. 



3 OBSERVATIONS 

As mentioned in the introduction, not all of the telescopes 
in the LBA network are capable of observing in the typi- 
cal geodetic/astrometric mode of dual S/X frequency bands 
with multiple spaced subbands and recorded to a Mark5 
VLBI system. 

The HOBART26, PARKES, and DSS45 have a MarkS 
recording system and a Mark4 (| Whitney et al.l [20041 ) base- 
band conversion rack. However fo r the first epoch the 
LBADR system (|Phillips et al.ll2009l ) was used. For all sub- 
sequent experiments the Mark5 system was used. Data were 
recorded onto normal Mark5 diskpacks and shipped to the 
Max-Plank Institut fiir Radioastronomie in Bonn for pro- 
cessing. 

The ATCA-104, MOPRA and CEDUNA only have the stan- 
dard LBA VLBI backend consisting of an Australia Tele- 
scope National Facility (ATNF) Data Acquisition System 
(DAS) with an LBADR recorder. The ATNF DAS only al- 
lows two simultaneous intermediate frequencies (IFs): either 
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Table 2. Observation Summary for each observing session 



Teleseope Recorder Polarization 



Frequency Bands 
(MHz) 
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RCP only 
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CEDUNA 
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SOI fi 
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8248- 


-8264 
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DSS45 




8456- 
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8472- 


-8488 






8488- 


-8504 


8504- 


-8520 


CEDUNA 


LBADR 


RCP only 


8200- 


-8216 


8216- 


-8232 






8456- 


-8472 


8472- 


-8488 


ATCA-104 


LBADR 


RCP & LCP 
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-8232 


MOPRA 
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-8472 
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-8488 
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-8488 
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-8536 


CEDUNA 


LBADR 
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8200- 


-8216 


8216- 


-8232 
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-8488 
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Experiment v271c, July 04 2009 
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2 frequencies or 2 polarizations. For each of these IFs the in- 
put 64 MHz analog IF is digitally filtered to 2 contiguous 
16 MHz bands. ATCA-104, MOPRA, and PARKES have two 
ATNF DAS, however the IF conversion system at each tele- 
scope means it is impractical to run in any modes other 
than 2 frequencies and dual polarization. The LBADR data 
format is not compatible with the Mark4 data processor at 
Bonn. A custom program was written to translate the data 
to MarkSB format then it was electronically copied using 
the Tsunami-UDP application to Bonn, before being copied 
onto Marks diskpacks. 

In the last experiment PARKES recorded on the LBADR 
system in parallel with the Mark5. As this included left cir- 
cular polarization (LCP) which was not recorded on the 
Marks system and since ATCA-104 and MOPRA recorded both 



RCP and LCP, the LBADR data were also sent to Bonn and 
LCP correlated against ATCA-104 and MOPRA. 

The main limiting factor for frequency selection was the 
LBADR backend at ATCA-104, MOPRA and CEDUNA where a 
setup with bands centered on ~8.2 and~8.5 GHz was cho- 
sen. For the telescopes with Mark5 recorders, the setup was 
chosen to overlap in frequency with the LBADR setup but 
also including more frequencies to gain sensitivity. The fre- 
quency setup was changed between experiments in order to 
explore the feasibility of recording at 512 Mbit/s at those 
stations that can support it. The setup for each observing 
experiment is described in Tabled 

The Australia Telescope Compact Array consists of six 
22 metre antennas and may observe as a single antenna or as 
a phased array of 5 dishe^. Since no tests of using a phased 
ATCA array as an element of the VLBI network for absolute 
astrometry observations were made before 2010, we decided 
to use a single dish of the ATCA in order to avoid the risk 
of introducing unknown systematic errors. 

The NASA Deep Space Network station DSS45 observed 
only 4.5 hours in v271a and 6.5 hours in v271b during inter- 
missions between receiving the telemetry from Mars orbiting 
spacecraft. 



4 CORRELATION 

The Bonn Mark4 Correlator was chosen for the correlation 
of these experiments for two primary reasons. Firstly, the 
correlator was extensively tested for use in space geodesy 
and absolute astrometry mode during 2000-2010 and is 
known to produce highly reliable group delays. Secondly, 
it was equipped with four Mark5B units and eight Mark5A 
units, which was a convenient combination since data from 
PARKES were recorded in the Mark5B format and data 
from ATCA-104, MOPRA, and CEDUNAoriginally recorded in 
LBADR format were transformed to the Mark5B format be- 
fore correlation. 

Preparation for correlation required about two months 
due to complications that arose from the scheduling of dif- 
ferent patching and channel outputs (fan out) for use at the 
stations. The stations that had Mark4 data acquisition racks 
delivered detailed log files, from which we could reconstruct 
the channels used, but the LBA stations did not provide log 
files and the track assignments had to be searched for by 
trial and error, which was a time-consuming process. The 
absence of log files also required some custom programming 
to reconstruct the log VEX file (Ivex), which the Mark IV 
correlator required to perform the correlation. Further, the 
setup that was chosen was incompatible with the capabil- 
ities of the hardware correlator and required an extensive 
work-around at the correlator. 

We chose to correlate with a window of 128 lags (cor- 
responding to a delay window width of 8 /xs) instead of the 
32 lags (2 ^s) normally used for stream correlation. This was 
to allow for potentially large clock offsets due to instrumen- 
tal errors, and for large a priori source coordinate errors. 
The integration time was chosen to be as small as possible 

^ The sixth dish is located at the fixed pad far away from other 
antennas, which makes its phasing with the rest of the compact 



array too difficult. 
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(0.5 s to 1.0 s) to allow for potentially large residual fringe 
rates due to source position errors. 

The Mark4 correlator can cope with only four stations 
simultaneously when correlating 128 lags with short inte- 
gration time. Since the experiment was observed with five 
to seven stations in right circular polarization (RCP) and 
two stations (atca-104 and mopra) also in left circular po- 
larization (LCP), the correlation had to be split in passes. 

As a by-product of the multi-pass correlation, some 
baselines were correlated more than once and this redun- 
dancy was exploited to select the correlation with the best 
SNR for each duplicated scan. For the LCP correlation, only 
one pass was required since only the LBA stations recorded 
LCP. PARKES recorded LCP only during v271c and this was 
enabled by the use of two backends (Mark4 and LBA) in par- 
allel. Those stations were re-correlated in a second pass to 
produce all four polarization products (left-left, right-right, 
left-right and right-left). For this pass we had to restrict the 
number of lags to 32 due to correlator hardware constraints, 
but the integration time was kept at 0.5 s. 

Fringe fitting was performed at the correlator using 
software program fourfit, the baseline-based fringe fit of- 
fered within the Haystack Observatory Package Software 
(hops) to estimate the residual delay upon which the post- 
correlation analysis was based. Inspection of fringe-fitted 
data is a convenient means of data quality control, as one 
can examine the correlated data on a scan-by-scan basis, 
looking for problems that might have occurred during cor- 
relation or at the stations and permits flagging of bad data 
immediately. 



5 DATA ANALYSIS 

The correlator either generates the spectrum of the cross- 
correlation function directly, or it can be easily derived 
from its raw output. The spectrum was processed with 
the software program fourfit, which for each scan and 
each baseline determined the phase delay rate, narrow-band 
delay and wide-band group delay (sometimes also called 
multi-band dela y ) tha t maximized the fringe amplitude. 
iTakahashi et all ((2OO0I ') give a detailed description of the 
fringe searching process and the distinction between narrow- 
band and wide-band group delays. The wide-band delay is 
more precise than narrow-band delay. The formal uncertain- 
ties of these delays, On and ct™ are computed by fourfit 
the following way: 

Vvl 

^ 27rA/ SNR 

" 27ro-/ SNR' 

where A/ the IF bandwidth, aj the dispersions of cyclic fre- 
quencies across the band, and SNR is the ratio of the fringe 
amplitude from the wide band fringe search to the rms of the 
thermal noise. The ratio of cr„/(T„ was in the range 27-28 
for the LCSl experiments, which gives for observations with 
typical SNR=30 ~ 40 ps and (t„ around one nanosecond. 

Of the 421 observed targets, 410 were detected. The list 
of 11 undetected sources is given in Table [S] In addition 111 
calibrators were observed, all of which were detected. 



Table 3. The list of target sources that were not detected in LCSl 
VLBI observations. The last three data columns contain the flux 
density in mjy extrapolated to 8.3 GHz (F), the spectral index 
(a), and comment (C). Value -9.99 means the spectral index was 
not available. Comments: 1) Not in AT20G catalogue; 2) Flagged 
as an LMC source in AT20G catalogue; 3) Galactic planetary 
nebula; 4) Flagged as an extended source in AT20G catalogue. 
The table is also available in the electronic attachment. 



J2000- 


name 


B1950 




RA 




Dec 
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J0404- 


-7109 


0404- 


-712 


04 
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00.99 


-71 


09 


09.7 


172 


-0.57 


1 


J0538- 


-6905 


0539- 


-691 


05 


38 


45.66 


-69 


05 


03.1 


204 


-9.99 
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J0552- 
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0551- 
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05 


52 


36.18 


-53 
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32.4 


182 


-0.34 
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38 


47.20 


-60 
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28.7 


196 


-0.08 


3 
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58 


02.92 


-57 


57 


42.6 


383 


0.21 


3 


JllOO- 


-6514 


1058- 


-649 


11 


00 


20.09 


-65 


14 


56.4 


179 


-0.08 


3 


J1150- 


-5710 


1147- 


-569 


11 


50 


17.87 


-57 


10 


56.0 


349 


0.23 


3 


J1325- 


-4302 


1322- 


-427 


13 


25 


07.35 


-43 


02 


01.8 


327 


-9.99 


4 


J1353- 


-6630 


1350- 


-662 


13 


53 
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5.1 Data analysis: source position determination 

The most challenging part of data analysis was resolving 
group delay ambiguities. The algorithm for fringe fitting im- 
plemented in FOURFIT searches for a global maximum in the 
Fourier transform of the cross-correlated function, averaged 
over individual IFs after correcting phases for a fringe de- 
lay rate and a narrow-band group delay. Fringe spectrum 
folding results in a rail of maxima of the Fourier transform 
of the cross-correlated function with exactly the same am- 
plitude and with spacings reciprocal to the minimum differ- 
ence in intermediate frequencies, 1/16.0 ■ 10® Hz — 62.5 ns. 
Within one half of the 62.5 ns range, the fringe spectrum has 
several strong secondary maxima. The second maximum at 
1/256.0-10** Hz ^ 3.9 ns has the amplitude 0.981 of the main 
maximum, the third maximum at 2/256.0 • 10* Hz ~ 7.8 ns 
has the amplitude 0.924. This high level of secondary max- 
imum amplitudes is due to our choice of intermediate fre- 
quencies that was determined by the hardware limitations. 
Due to the presence of noise in the cross-correlated func- 
tion, the difference between the fringe amplitude at the 
main maximum and at the secondary maxima is random, 
and therefore, the group delay is determined with an uncer- 
tainty ±1 — 3 of the spacing between the secondary maxima, 
3.9 ns. 

It should be noted that typical group delay ambiguity 
spacings in geodetic observations are in the range 50-200 ns. 
For successfully resolving ambiguities, as a rule of thumb the 
predicted delay should be known with an accuracy better 
than 1/6 of the ambiguity spacing, 600 ps in our case. This 
is a challenge. The random errors of the narrow-band group 
delays are too high to use directly for resolving ambiguities 
in wide-band group delays (see Fig. [2]). In the framework 
of traditional data analysis of geodetic VLBI observations, 
the maximum uncertainty in prediction of the path delay 
is due to a lack of an adequate model for the path delay 
in the neutral atmosphere. These errors are in the range 
of 300-2000 ps. The path delay through the ionosphere at 
8.3 GHz is in the range 30-1000 ps. In addition, an error 1" 
in a priori source position (a typical AT20G position error) 
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causes an error in a priori time delay of up to 30 ns on a 
1700 km baseline. 



5.1.1 Group delay ambiguity resolution 

However, it is premature to conclude that resolving group 
delay ambiguities is impossible. Firstly, we need to use 
a state-of-the art a priori model. Our computation of 
th eoretical time delays in genera l follow s the approach 
of i overs. Fanselow fc Jacobs C.Sl l|l998l ) with some re- 
finements. The most significant ones are the foUow- 
ing. The advanced expre ssion for time delay derived by 
iKopeikin fc Schafeil (jl999f l in the framework of general rel- 
ativity was used. The displacements caused by the Earth's 
tides were computed using the num erical values of t he gen- 
eralized Love numbers presented bv lMathewsl (1200 ll) follow- 
ing a rigorous algorithm described by Petrov fc Mai (|2003l ) 
with a truncation at a level of 0.05 mm. The displacements 
caused by ocean loading were computed by convolving the 
Greens' functions with ocean tide models. Th e GOT99. 2 
model of diurnal and semi-diurnal ocean tides (lRavlll999t) 



10 



Ooj I2OOOI') 



the NA099 model jMatsumoto. Takanezawa" 

of oc ean zonal tides, the equilibrium model ( Petrov fc Mai 
|2003| ) of the pole tide, and the tide with period of 18.6 
years were used. Station displacements caused by the at- 
mospheric pressure loading were computed by convolving 
the Greens' f unctions that describe the elastic properties 
of the Earth (|Farrelll Il972l ) with the output of the atmo- 
spher e NCEP Reanalysis numerical model (|Kalnav et al.l 
1 19961 '). T he algorithm o f computations is described in full 
detail in IPetrov fc Bovl (|2004 ). The displacements due to 
loading caused by variations of soil moist ure and snow cover 
in accordance with GLDAS Noah model (|Rodell et al.ll2004l ') 
with a resolution 0.25° x 0.25° were computed using the same 
technique as the atmospheric pressure loading. The empir- 
ical model of harmonic variations in the Earth orientation 
parameters heo_20101111 derived from VLBl observa tions 
according to the method proposed by IPetrovl (|2007l ) was 
used. The time series of UTl and polar motion derived by 
the NASA Goddard Space Flight Center operational VLBI 
solutions were used a priori. 

The a priori path delays in the neutral atmosphere 
in the direction of observed sources were computed by 
numerical integration of differential equations of wave 
propagation through the heterogeneous media. The four- 
dimensional field of the refractivity index distribution was 
computed using the atmospheric pressure, air temperature 
and specific humidity taken from the output of the Modern 
Era Retros pective- Ana lysis for Research and Applications 
(MERRA) (jSchubert et al. 2008). That model presents the 
atmospheric parameters at a grid 1/2° x 2/3° x 6'' at 72 
pressure levels. 

Secondly, we made a least square (LSQ) solution using 
narrow-band delay. The positions of all target sources were 
estimated, as well as clock functions for all stations except 
the one taken as a reference and the residual atmosphere 
path delay in zenith directions. The clock function was mod- 
eled as a sum of the 2nd degree polynomial over the exper- 
iment and B-spline of the first order with the time span 
60 minutes. The residual atmosphere path delay in zenith 
direction was modeled with B-spline of the first order with 
time span 60 minutes. Constraints on rate of change of clock 
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Figure 2. The post-fit residuals of narrow-band group delay LSQ 
solution at the baseline CEDUNa/hobart26 in experiment v271a 
on February 05, 2008. The bar in the middle of the plot corre- 
sponds to the ambiguity spacing of the wide-band path group 
delay. 
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Figure 3. The wide-band a priori group delays after substitution 
of adjustments to clock function, atmosphere path delay in zenith 
direction and source positions from the narrow-band LSQ solution 
at the baseline CEDUNA/parkes in experiment v271a on February 
05, 2008. The bar in the left side of the plot corresponds to the 
ambiguity spacing of the wide-band path group delay. 



function and atmosphere path delay were imposed. After 
the removal of outliers (1-2% of points), the weighted root 
mean square (wrms) of residuals was 1-4 ns. An example of 
narrow-band postfit residuals is shown in Fig. [2| 

The next step was to substitute the adjustments to 
clock function, atmosphere path delay in zenith direction 
and source positions from the narrow band LSQ solution 
to the differences between the wide-band path delays and 
the theoretical group delays. An example of enhanced wide- 
band group delays after the substitution and group delay 
ambiguity resolution is shown in Fig. [3l 

It is still difficult to resolve ambiguities at long base- 
lines, but relatively easy to resolve the ambiguities at the 
inner part of the array: ATCA-104, DSS45, mopra, parkes. 
The group delay ambiguity resolution process starts from 
the baseline with the least scatter of a priori wide-band de- 
lays. After group delay ambiguity resolution at the first base- 
line and temporary suppression of observations with ques- 
tionable ambiguities, the LSQ solution with mixed delays 
is made: wide-band group delays at baselines with resolved 
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Figure 4. The wide-band a priori group delays at the baseline 
CEDUNA/parkes in experiment v271a on February 05, 2008 after 
substitution of adjustments to clock function, atmosphere path 
delay in zenith direction from the narrow-band LSQ solution and 
source positions from the wide-band LSQ solution at other base- 
lines and group delay ambiguity resolution. The bar in the central 
side of the plot corresponds to the ambiguity spacing of the wide- 
band path group delay 



ambiguities and narrow-band group delays at other base- 
lines. In addition to other parameters, baseline-dependent 
clock misclosures are estimated. 

In the absence of instrumental delays, the differences 
between wide-band and narrow-band delays would be the 
zero mean Gaussian random noise. Instrumental delays in 
the analogue electronics cause systematic changes of these 
differences in time. Fortunately, these systematic changes 
are smooth and small enough to allow solving for ambiguities 
with spacings as small as 3.9 ns. Since the scatter of postfit 
residual of wide-band group delays is a factor of 10-50 less 
than the ambiguity spacings, when the instrumental delay 
is determined, the ambiguities can be resolved. Therefore, 
the feasibility of resolving wide-band group delays hinges 
upon the accuracy of determination of the differences be- 
tween narrow-band and wide-band group delays from anal- 
ysis of narrow-band delays. 

An example of the a priori wide-band group delay with 
the substituted adjustments to sources positions found from 
the wide-band group delay solution at other baselines is 
shown in Fig. |4l After resolving ambiguities for the most of 
observations, the points that were temporarily suppressed 
were restored and used in the solution. 



5.1.2 Ionosphere path delay contribution 

Single band VLBI data are affected by a variable path delay 
through the ionosphere. We attempted to model this path 
delay using GPS TEC maps provided by the CODE analy- 
sis ce nter for proce ssing Global Navigation Satellite System 
data (|Schaeil Il998l ). Details of computing the VLBI iono- 
spheric path delay using TEC maps from GPS are given in 
iPetrov et al.l (|201ll ). The model usually recovers over 80% 
of the path delay at baselines longer than several thousand 
kilometres. However, applying the ionosphere TEC model 
to processing LCSl experiment did not improve the solution 
and even degraded the fit. The same model and software pro- 
gram certainly improved the fit and improved results when 
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Figure 5. The difference in rms of postfit residuals with respect 
to the reference dual band X/S solution for two trial solutions: 
X-band only data (solid discs) and X-band data with ionosphere 
path delay from GPS TEC map applied (circles). 



applied to processing observations on intercontinental base- 
lines. 

In order to investigate the applicability of the reduction 
for the path delay in the ionosphere based on GPS TEC 
maps, we processed 29 IVS dual-band geodetic experiments 
that included the 1089 km long baseline hobart26/parkes. 
We ran three solutions that included the data only at this 
baseline. The first reference solution used ionosphere-free 
linear combinations of X-band and S-band observables. The 
second solution used X-band only group delays. The third 
solution used X-band only data and applied the reduction 
for the path delay in the ionosphere from GPS TEC maps. 

We discarded two experiments: one had a clock break 
at HOBART26 and another had the rms of postfit residual 
a factor of 5 greater than usually due to a warm receiver. 
The rms differences of postfit residuals with respect to the 
reference dual band X/S solution for two trial solutions are 
shown in Fig. [S] one with X-band only data and another 
with X-band data with ionosphere path delays from GPS 
TEC map applied. Considering the reference solution based 
on ionosphere-free group delays as the ground truth, we ex- 
pected that the rms differences between the X-band only 
solution (denoted with circles in Fig. [5]) be positive, and the 
rms differences of the X-band solution with the ionosphere 
path delay from GPS TEC maps applied be also positive 
but less than the rms differences between the X-band only 
solution. Instead, we see that the rms of postfit residuals of 
the X-band only solution are less than the rms of dual-band 
reference solution and applying the ionosphere path delay 
from GPS TEC maps does not affect the rms significantly. 

We then computed the rms of the contribution of the 
ionosphere path delay from dual-band X/S observations and 
the rms of the differences between the path delay in the iono- 
sphere from the dual-band X/S observations and the GPS 
TEC ionosphere maps. They are shown in Fig. (6) We see 
that in 2002-2004 the ionosphere path delay from GPS TEC 
maps was coherent with the dual-band ionosphere path de- 
lay estimate. In 2005-2010 the ionosphere path delay was 
small with rms ~30 ps and the ionosphere path delay from 
the GPS TEC model was only partially coherent with the 
ionosphere path delay estimates from dual-band VLBI ob- 
servations. The baseline length repeatability, defined as the 
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Figure 6. The rms of the ionosphere contribution from dual-band 
X/S observations (solid disks) and the rms of the differences of the 
path delay in the ionosphere from the dual-band X/S observations 
at the baseline hobart26/parkes and the GPS TEC ionosphere 
maps (circles). 
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Figure 7. The rms of the ionosphere contribution from dual-band 
X/S observations (solid disks) and the rms of the differences of the 
path delay in the ionosphere from the dual-band X/S observations 
at the 1088 km long baseline la-vlba/ov-vlba and the GPS 
TEC ionosphere maps (circles). 



rms of baseline length estimates after subtracting the secular 
drift due to tectonics, is the minimum when ionosphere-free 
linear combinations of X/S observables are used: 17.0 mm, 
grows to 17.4 mm when X-band group delay are used, and 
it is the maximum, 18.4 mm, when the ionosphere path de- 
lay reduction based on GPS TEC model is applied to the 
X-band observables. 

We compared these results with analogous observations 
in the northern hemisphere. We picked the baseline LA- 
vlba/ov-vlba that has almost exactly the same length as 
the baseline hobart26/parkes: 1088 km against 1089 km. 
We processed 88 experiments at this baseline in a similar 
fashion as we processed the hobart26/parkes baseline. Re- 
sults are shown in Fig.[7l Applying the ionosphere path delay 
from GPS has reduced the rms of post-fit residuals at this 
baseline. 

Our interpretation of this phenomena is that the iono- 
sphere path delay from GPS TEC maps at baselines 1-2 
thousands kilometres has the floor around 30 ps. The domi- 
nant constituent in the ionospheric path delay at these scales 
during the solar minimum are short-periodic scintillations 



that the GPS TEC model with time resolution 2 hours and 
spatial resolut ion 500 km does not adequate ly represent. It 
is also known (|Hernandez-Paiares et al"]|2009l ) that the accu- 
racy of the GPS TEC model is worse in the southern hemi- 
sphere than in the northern hemisphere due to the disparity 
in the GPS station distribution. 

Applying the GPS TEC map ionosphere path delay re- 
duction effectively added noise in the data. Therefore, we did 
not apply the ionosphere path delay in our final solution. 



5.1.3 Re-weighting observations 

According to the Gauss-Markov theorem, the estimate of 
parameters has the minimum dispersion when observation 
weights are chosen reciprocal to the variance of errors. The 
group delays used in the analysis have errors due to the 
thermal noise in fringe phases and due to mismodeling the 
propagation delay: 



2 2 2 2 



(2) 



where ath is the variance of the thermal noise, a„a, and <Jio 
axe the variances of errors of modeling the path delay in the 
neutral atmosphere and the ionosphere respectively. 

A rigorous analysis of the errors of modeling the path 
delay in the neutral atmosphere is beyond the scope of this 
paper. Assuming the dominant error source of the a priori 
model is the high frequency fiuctuations of the water vapor 
at time scales less than 3-5 hours, we sought a regression 
model for the dependence of a^a on the non-hydrostatic 
component of the slanted path delay through the atmo- 
sphere. We made several trial runs using all observing ses- 
sions under geodesy and absolute astrometry programs for 
30 years and added in quadrature to the a priori uncertain- 
ties of group delay an additive correction: 



2 _ 2 
""used — <^th 



(3) 



where is the contribution of the non-hydrostatic con- 
stituent of the slanted path delay, is the non-hydrostatic 
path delay in zenith direction computed by a direct inte- 
gration of equations of wave propagation through the atmo- 
sphere using the refractivity computed using the MERRA 
model, and a is a coefficient. We found that the baseline 
length repeatability defined as the rms of the deviation of 
baseline length with respect to the linear time evolution 
reaches the minimum when a = 0.02. Therefore, we adopted 
this value in our analysis of LCSl experiments. For typical 
values of t^, the added noise is 8 ps in zenith direction and 
16 ps at the elevation of 30°. 

We computed the rms of ionosphere contributions from 
dual band VLBI group delays and from GPS TEC maps. 
Fig. [8] shows the dependence of the square of the rms from 
dual-band VLBI versus the square of the rms from GPS. 
The slope of the regression straight line is 0.992±0.004. This 
dependence suggests that although the ionosphere path de- 
lay from TEC maps from GPS analysis is too noisy to be 
applied for reduction of observations at short baselines dur- 
ing the Solar minimum, it correctly predicts the variance of 
the ionosphere path delay. We assigned the variance of the 
mismodeled ionosphere path delay to the ionosphere contri- 
bution to delay computed from TEC maps: cjio = Tio. 
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Figure 8. The dependence of the ionosphere path delay cr^ from 
dual-band VLBI versus cr^ from TEC maps from GPS analysis. 
Data in time range 2005.0-2010.6 (Solar minimum) were used. 

We also computed, for each experiment and for each 
baseline, ad hoc variances of observables that after being 
added in quadrature make the ratio of the weighted sum of 
squares of post-fit residuals to their mathematical expecta- 
ti on close to un i ty. Th is computation technique is presented 
m IPetrov et al] (|201lf ). The ad hoc variance was applied to 
further inflate the formal uncertainties of the observables 
that have already been corrected for the inaccuracy of the 
a priori model of wave propagation through the ionosphere 
and neutral atmosphere (expression [2J . In contrast to cr^^ 
and (j^o, the baseline-dependent ad hoc variance is elevation 
independent. 

5.1.4 Global LSQ solution 

We ran a single global LSQ solution using all available dual- 
band VLBI observations under geodesy and absolute astron- 
omy programs from 1980.04.12 through 2010.08.04, in total 
7.56 million observations, and the X-band VLBI data from 
4 LCSl observing sessions. The RCP and LCP data were 
treated as independent experiments. The following param- 
eters were estimated over the global dataset: coordinates 
of 4924 sources, including 410 detected target objects in 
the LCSl campaign (see Fig. (9)1; positions and velocities 
of all stations; coefficients of expansion over the B-spline 
basis of non-linear motions for 17 stations; coefficients of 
harmonic site position variations of 48 stations at four fre- 
quencies: annual, semi-annual, diurnal, semi-diurnal; and 
axis offsets for 67 stations. In addition, the following pa- 
rameters were estimated for each experiment independently: 
station-dependent clock functions modeled by second order 
polynomials, baseline-dependent clock offsets, the pole co- 
ordinates, UTl-TAI, and daily nutation offset angles. The 
list of estimated parameters also contained more than 1 mil- 
lion nuisance parameters: coefficients of linear splines that 
model atmospheric path delay (20 minutes segment) and 
clock function (60 minutes segment). 

The rate of change for the atmospheric path delay and 
clock function between adjacent segments was constrained 
to zero with weights reciprocal to 1.1 ■ 10~^^ and 2 • 10^^** 
respectively in order to stabilize our solution. We apply no- 
net rotation constraints on position o f 212 sources m arked 
as "defining" in the ICRF catalogue (|Ma et all [19981 ) that 
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Figure 9. The distribution of sources with milliarcsecond posi- 
tions at the southern sky. Circles denote sources known by 2008. 
Filled green discs denote sources from the LCSl catalogue. The 
Galactic plane shown by red line. 



requires the positions of these source in the new catalogue 
to have no rotation with respect to the position in the ICRF 
catalogue. 

The global solution sets the orientation of the array 
with respect to an ensemble of ~ 5000 extragalactic remote 
radio sources. The orientation of that ensemble is defined 
by the series of the Earth orientation parameters evaluated 
together with source coordinates. Common 111 sources ob- 
served in LCSl as atmosphere and amplitude calibrators 
provide strong connection between the new catalogue to the 
old catalogue of compact sources. 



5.2 Error analysis of the LCSl catalogue 

To assess the systematic errors in our results we exploited 
the fact that 111 known sources were observed as amplitude 
and atmospheric calibrators. Positions of these sources were 
determined from previous dual-band S/X observations with 
accuracies better than 0.2 mas. We sorted the set of 111 cal- 
ibrators according to their right ascensions and split them 
into two subsets of 55 and 56 objects, even and odd. We ran 
two additional solutions. In the first solution we suppressed 
55 calibrators in all experiments but LCSl and determined 
their positions solely from the LCSl experiment. In the sec- 
ond solution we did the same with the second subset. Con- 
sidering that the positions of calibrators from numerous S/X 
observations represent the ground truth, we treated the dif- 
ferences as LCSl errors. 

We computed the x^/ndf statistics for the differences 
in right ascensions and declinations Ac and and sought 
additional variances Va and vs that, being added in quadra- 
ture to the formal source position uncertainties, aa,i and oj^i. 
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make them close to unity: 

Aa^ cos^ Si 



X5 



(4) 



We found the following additive corrections of the un- 
certainties in right ascensions scaled by cos 5 and for decli- 
nations respectively: Vc, — 1.44 mas and vs=0.51 mas. We 
do not have an explanation why the additive correction for 
scaled right ascensions is 3 times greater than for declina- 
tions. After applying the additive corrections, the wrms of 
source position differences are 1.8 mas for right ascension 
scaled by cos 5 and 1.5 mas for declinations. 

The final infiated uncertainties of source positions, 



;(/) and a^if), are 



0"q + l-'a I COS^ 5 



(5) 



5.3 Data analysis: correlated flux density 
determination 

Each detected source has from 3 to 60 observations, with 
a median value of 25. Imaging a source with 25 points at 
the uv plane is a difficult problem, and the dynamic range 
of such images will be between 1:10 and 1:100, which is 
far from spectacular. Images produced with the hybrid self- 
calibration method will be presented in a separate paper. In 
this study we limited our analysis to mean correlated fiux 
density estimates in three ranges of lengths of the baseline 
projections onto the plane tangential to the source, without 
inversion of calibrated visibility data. Information about the 
correlated flux density is needed for evaluation of the re- 
quired integration time when an object is used as a phase 
calibrator. 

First, we calibrated raw visibilities v for the a priori 
system temperature Tsys and antenna gain G(e): Acorr ~ 
V ■ Tsys{t, e) /G{e). The coefficients of antenna gain expan- 
sions into polynomials over elevation angle e are presented in 
Table |4l System temperature was measured at each station, 
each scan. 

At the second step, we adjusted antenna gains us- 
ing publicly available brightness distributions of calibrator 
sources made with observations under other programs. We 
compiled The Database of Brightness Distributions, Corre- 
lated Flux Densities and Images of Compact Radio Sources 
Produced with VLSfl from authors who agreed to make 
their imaging results publicly available. Among 111 sources 
observed as calibrators, images of 14-27 objects at each ex- 
periment were available. These a re images in the f orm o f 
CLEAN components mainly from lOiha et all (|2004 l2005h . 
CLEAN components of source brigh tness distributions from 
analysis of the TAN AMI program (|Oiha et all I2OIOI ') that 
observed with the LB A concurrently with the LCSl were 
not available, but the parameters of one-component Gaus- 
sian models that fit the core regions were published. For 

Available at | http : //astrogeo ■ org/vlbi_lmages| 



those sources for which both a set of CLEAN components 
and the parameter of the Gaussian one-component model 
were available, we used CLEAN components. 

We predicted the correlated flux density Fcorr for each 
observation of an amplitude calibrator with known bright- 
ness distribution, as 



(" x + v y) 



(6) 



± corr 



5" g 4ln 2 ^^"^ ''^^ 6-\~v sin d)"^ '^b'^ ( — u sin Q-^-v cos Q)"^) 



where d is the correlated flux density of the ith CLEAN 
component with coordinates x and y with respect to the 
center of the image; u and v are the projections of the base- 
line vectors on the tangential plane of the source; and a and 
6 are the FWHM of the Gaussian that approximates the 
core, and 9 is the position angle of the semi-major axis of 
the Gaussian model. 

Then we built a system of least square equations for 
all observations of calibrators with known images used in 
astrometric solutions: 



p — 

-t corr — 



(7) 



and after taking logarithms from left and right hand sides 
solved for corrections to gains g for all stations. Finally, 
we applied corrections to gain for observations of all other 
sources. 

The correlated fiux density is a constant during observ- 
ing session only for unresolved sources. For resolved sources 
the correlated fiux density depends on the projection of the 
baseline vector on the source plane and on its orientation. 
We binned the correlated fiux densities in three ranges of the 
baseline vector projection lengths, 0-6 MA, 6-25 MA, 25-50 
MA, and found the median value within each bin. These bins 
correspond to scales of the detected emission at > 30 mas, 
7-30 mas, and < 7 mas respectively. 

The uncertainties of our estimates of the correlated 
density depend on the thermal noise described as ath = 
Jcorr/SNR and errors of gain calibration. The variance of 
the thermal noise was in the range of 1 to 6 mjy depend- 
ing on the sensitivity of a baseline, with the median value 
of 3 mJy. The LSQ solution for gains provided the variance 
of the logarithm of gains. Assuming the calibration errors 
to be multiplicative in the form g{l + e), where g is gain 
and e is the Gaussian random variable, we can evaluate the 
contribution of gain errors in the multiplicative uncertainty 
of Fcorr + fJ'), where fi the Gaussian random variable. Its 
variance is evaluated as 



cr(M) 



exp 



^i-^Cov(si, Si) + Cov(s2, 52) + 2Cov(si, S2)^ 



(8) 



where Cov(si, Sj) is the covariance of the logarithm of gain 
between the i-th and j-th station. Multiplicative gain uncer- 
tainties are in the range 0.08-0.1 for DSS45 and 0.02-0.05 for 
other stations. The gain uncertainties for DSS45 are higher, 
because it observed 3-5 times less than other stations. The 
total variance is a sum in quadrature of Uc and ath- We 
should note that our estimate of systematic errors does not 
account for possible errors in gain curve determination. A 
systematic error in the gain curve would directly affect our 
estimate of the correlated fiux density. It will also affect the 
maps of calibrators sources that we took from literature, and 
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k=5 

Table 4. The coefficients of antenna gain polynomials over elevation ang le expressed in radians. Gain = DPFU ^ a^.E''. The DFPU 

k=o 

is shown in the second columns. Columns 3-8 hold coefficients of degree through 5. 



Station DFPU K/Jy The coefficients of polynomial for gain as a function of the elevation angle in degrees 
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Figure 10. The probability density histogram of correlated flux 
densities in LCSl experiments. The dotted line with dicks denotes 
correlated flux densities at baseline projection lengths 0-6 MA, 
the solid line with discs denotes the correlated flux density at 
projection lengths 25-50 MA, 



thus indirectly affect our estimates of gain corrections. We 
do not have information about uncertainties of gain curves 
of LB A antennas. 

Table [5] displays 12 out of 633 rows of the catalogue of 
correlated flux density estimates. The full table is available 
in the electronic attachment. Column 1 and 2 contain the 
J2000 and B1950 lAU names of a source; column 3 contains 
the observation date of the experiment; the 4th column con- 
tains the status of the source: Cal if it was used as amplitude 
calibrator; the 5th columns contains the number of RCP ob- 
servations used in processing. Columns 6, 7, and 8 contain 
median values of correlated flux densities determined in that 
experiment at baseline projection lengths 0-6 MX, 6-25 MX, 
and 25-50 MA respectively. Columns 9, 10, and 11 contain 
the arithmetic mean of correlated flux density uncertainties 
that accounts for both thermal noise and uncertainties in 
calibration. If no data were collected to that range of base- 
line projection lengths, —1.0 is used as a substitute. Column 
12 contains the experiment name. 

Fig. [To] shows the probability density histogram of cor- 
related flux densities in LCSl experiments. 
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Figure 11. The histogram of the semi-major axes of inflated 
position error ellipses among 410 sources in the LCSl catalogue. 



6 THE LCSl CATALOGUE 

Table |6] displays 12 out of 410 rows of the LCSl catalogue of 
source positions. The full table is available in the electronic 
attachment. Column 1 and 2 contain the J2000 and B1950 
lAU names of a source; column 3, 4, and 5 contain hours, 
minutes and seconds of the right ascension; columns 6, 7, 
8 contain degrees, minutes and arcseconds of declination. 
Columns 9 and 10 contain inflated position uncertainties in 
right ascension (without multiplier cos 5) and declination in 
milliarcseconds. Column 11 lists the correlation coefficient 
between right ascension and declination, and column 12 con- 
tains the total number of observations used in position data 
analysis, including RCP and LCP data. Columns 13, 14, 
and 15 contain the median estimates of the correlated flux 
density in Jansky over all experiment of the source in three 
ranges of baseline projection lengths: 0-6 MX, 6-25 MX, 
and 25-50 MX. Columns 16, 17, and 18 contain estimates 
of the correlated flux density uncertainties in Jansky. 

Of 421 sources observed, 410 were detected in three or 
more scans. In four LCSl observing sessions, 17,731 obser- 
vations out of 19,494 were used in the LSQ solution together 
with 7.56 million other VLBI observations. The semi-major 
error ellipse of inflated position uncertainties varies in the 
range 1.4 to 16.8 mas with the median value of 2.6 mas. The 
distribution of sources on the sky is presented in Fig. [9] The 
histogram of position errors is shown in Fig. 111! 
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Table 5. The first 12 rows of the catalogue of source correlated flux density estimates for both 410 target sources and 111 calibrator 
sources in each experiment. Some sources were observed on more than one experiment. The table columns are explained in the text. The 
full table is available in the electronic attachment. 
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Table 6. The first 12 rows of the LCSl catalogue of source positions of 410 target sources. The table columns are explained in the text. 
The full table is available in the electronic attachment. 



J2000-name B1950-namc Right ascension Declination ""(cf) o'{5) Corr #Obs Correl. Flux Correl. Flux unccr. 
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7 SUMMARY 

The absolute astrometry LBA observations turned out 
highly successful. The overall detection rate was 97% — the 
highest rate ever achieved in a VLBI survey. If we exclude 
extended sources, non-AT20G sources and the six planetary 
nebulae included in the candidate list due to an oversight, 
the detection rate is 99.8%! We attribute this high detection 
rate to two factors. Firstly, the AT20G catalogue is highly 
reliable and is biased towards very compact objects. Select- 
ing candidates based on the simultaneous AT20G spectral 
index proves to be a good methodology. Secondly, the LBA 
has very high sensitivity. The baseline detection limit over 
2 minute of integration time varied from 7 mjy to 30 mjy, 
with 7-12 mJy at baselines with parkes. 

We have successfully resolved group delay ambiguities 
with spacing 3.9 ns for all observations. This became possible 
using the innovative algorithm exploiting relatively low level 
of instrumental group delay errors. 

We have determined positions of 410 target sources 
never before observed using VLBI, with a median uncer- 
tainty of 2.6 mas. Error analysis showed a moderate contri- 
bution of the mismodeled ionosphere path delay to the over- 
all error budget. Both random and systematic errors are ac- 
counted for in the uncertainties ascribed to source positions 



exploiting the overlap of 111 additional sources observed in 
LCSl experiments with their positions known from prior ob- 
servations. The positional accuracy of the LCSl catalogue 
is a factor of 350 greater than the positional accuracy of the 
AT20G catalogue, corresponding to the ratio of the maxi- 
mum baseline lengths of the LBA and the ATCA. The new 
catalogue has increased the number of sources with declina- 
tions < —40° from 158 to 568, i.e. by a factor of 3.5. 

We determined correlated flux densities for 410 target 
and 111 calibrator sources, and presented their median val- 
ues in three ranges of baseline projection lengths. The cor- 
related flux density of the target sources varied from 0.02 
to 2.5 Jy, and was in the range 80-300 mJy for 70% of the 
sources. The uncertainties of the correlated flux densities are 
estimated to be typically 5-8%. 

This observing program is continuing. By November 
2010, four additional twenty four hour experiments had been 
observed with several more observing sessions planned. 
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